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On the Propagation of an Arbitrary Electro -Magnetic 

Disturbance, on Spherical Waves of Light and 

the Dynamical Theory of Diffraction. 

By Professor H. A. Kowland. 



INTRODUCTION. 

In the year 1849 the great paper of Stokes "On the Dynamical Theory of 
Diffraction " was read before the Cambridge Philosophical Society, and this has 
remained until the present day the standard upon this important subject. 

The method of Stokes was based upon the old elastic solid theory of light, 
and gave the following conclusions : 

First. That when the incident light was plane polarized, the diffracted light 
from a small orifice was also plane polarized in such a manner that the displace- 
ment was in the same plane as that of the medium at the orifice. So that if a 
sphere was drawn with the orifice as a center and meridians drawn on the sphere 
with the axis in the direction of the vibration at the orifice, then these meridians 
represented the direction of displacement in the diffracted light. 

Second. The intensity of the polarized light was represented as follows : 
Let h represent the angle between the incident ray prolonged and the diffracted 
ray, and let <p be the angle between the diffracted ray and the direction of the 
displacement at the orifice. Then the intensity of the diffracted light around a 
very small orifice will be proportional to 

(1 + cos &Y sin 3 ^). 
The presence of the term in <£> indicates that the intensity of the diffracted light 
at a given point varies as one rotates the plane of polarization. In ordinary 
light the intensity varies as 1 + cos 3 S . 

Stokes and others have attempted to determine the relation between the 
plane of polarization and the direction of displacement by means of the first 
relation, but they have not agreed with one another, and this want of agreement 
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has usually been assigned to the fact that the gratings used have been ruled on 
glass rather than in free space, as the equations indicate. 

On examining this question from the point of view of the electro-magnetic 
theory of light, I have been led to entirely different results. But as the elastic 
solid theory for an incompressible solid must agree with the electro-magnetic 
theory, I have been led to examine the theory of Stokes, and believe that I have 
now discovered an error, which, if it were corrected, would lead to my result. 

The results which I reach are as follows : 

First. The plane of polarization of the diffracted light is determined as 
follows : Draw a sphere around the orifice and mark the point on the sphere 
where the incident light enters it. Through this point draw circles on the sphere 
whose planes are parallel to the electrostatic displacement at the orifice, and 
these. circles give the direction of the electrostatic disturbance in the diffracted 
beam. Repeat the same for the magnetic disturbance, which is at right angles 
to the electrical disturbance, and the circles indicate the direction of the magnetic 
disturbance in the diffracted beam. These two systems of circles are orthogonal 
to one another, as they should be. 

Second. The intensity of the diffracted light around a very small orifice is 
symmetrical around the incident ray prolonged and is proportional to 

(1 + cos of, 

and the same expression applies to ordinary light. 

Although the theory of diffraction forms the most interesting part of my 
paper, yet I have thought it worth while to treat of the general problem of 
spherical waves of light, which I have not seen considered anywhere else. The 
method is similar to that used in sound and the theory of heat conduction. 

Sphebioal Waves. 

Let F , G , H ; F u G x , etc., be a system of vectors derived from one 
another by the equations 

m+1— dy dz 

(rro + 1- ~dz~~~dV 

TT — dGm dFm 
m+1 ~ dx dy ' 
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Then, as is well known, all this system of vectors necessarily satisfy the equation 

of continuity 

dF m . dO m dH m 

dx dy dz ' 

except the primitive ones F Q , G , H , and even these can be made to satisfy the 

equation by the addition of terms of the form -7- , — , — • 

\Xn<c ^y *** 

The equations of light-waves either on the electro-magnetic theory or the 
elastic solid theory are of the form 

dtr ( dz ) dz 

dF dO dH 
dx dy dz 

In the electro-magnetic theory v is always zero. In the elastic solid theory 
the terms in v give the wave of normal disturbance like that of sound, and can 
nearly always be omitted, as the transverse and normal waves are immediately 
separated when v is different from V, as it always must be. As we wish to use 
the derived functions only, the terms in J" can also be omitted, as Maxwell has 
shown, and hence we can write on either theory 

<PF 



dP 
d?G 



= V*A 2 F 



(PTT 

mr 

V indicates the velocity of a plane wave in the medium. 

Bach one of these equations is the same as that of sound and can be solved 

in a similar manner. If we have solutions of these equations F , G , H , which 

do not satisfy the equation of continuity, then we can get other solutions 

F m , G m , H m from these which shall satisfy this condition. Denote a/ — 1 by i 

and let V_ in+1) be a solid harmonic. Then as F_ (n+:) is a homogeneous function 
Vol. vi. 
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of x, y and z of the degree — (n + 1), we have 



x 



dV-(n + l) ■ ^^-(tt + l) | ^rff-^-fl) 



+y ^^ ±i ) +2 ^^ ±1 ) == _ (w+1)7 _ (n+1) 



da; u dy dz 

and A 2 F_ (n+1) = 0. 

Make i^ = C , n p*F_ (n+1)e <°~ lw ' ) ~ F0 > 

where (7 is a function of p of a complex form 

C n = A n —iB n 

p 2 =x 2 + 2/ 2 +2 8 - 
This form has the advantage that by the addition of another term of the same 
form in which i has a negative value, the sum will reduce to a real form with 
circular arcs instead of powers of s ; a is introduced for generality. Substituting 
this form we have for the determination of G,* the equation 



cpa 



* + 2(a-^-^<7„=0. 



dp* v 'dp p 

Writing c-=a — ib we readily find, as others have found before for this 
equation, 

f n n + 1 n(n»-l») »+2 rcfa 8 - l>*-2») ] 

C n - C jl - y— + 2.4 ~cy 2X6 + et<5 ' J 

The series ends after n + 1 terms, and so we can put it in the form 
^=(-ir(l.3.5...(2n-l))^{l- C p + A^ c2p8 

_ _2,2_ (n-lX»-2) 33 | 

1.2.3 (2n-l)(2»— 2) C P + etc -J 

The following properties of the functions C n are useful. "Writing c = a — ib, 
we have n _ 9 <, n n x 

Making C„ = X»„p i £ -p(o-<i,) , 

we obtain cPD. , 1 dfl„ /, .... , (w+j)j \_ n 

If we replace p by a new variable equal to i (a — ib) p , this reduces to the equation of Bessel's func- 
tions : hence we can write D„ — J, n ■ ^(6 + ia) p , 

with similar terms for Go and flo . 

Hence C„ = C.V* -p(o ~ <B) .7' (B _ M) (& + «i)p, 

where CV is a function of w to be determined, and one must add a corresponding term with — i in 

place of + i . 
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ep / c dp 



o^=o.(i- n -^)+ ldc - 



e"' 



Putting cp = cpVl + 2s we have 

Extending the notation so that we should write C n (p) in the place of C n as above, 
we have, if we put pes = q 

^ e '=f^ { °.(p) + f °^p) + o a -+»^) + etc -} 

These expansions are useful in many calculations from the equations, and I 
believe they are given here for the first time, together with the differential 
expression from which they are derived. 

Hence we have obtained a complete solution of one of the equations. Now 

the general value of V_ in+1) is 

d* 1 

K.(» + D= (— l Y ^ dht,... dh n ~J' 

By varying the axes \, h^, h 3 , etc., the value and form of V_ {n+1) will change, 

and can thus have many values while G remains unaltered. Select three of 

these forms and- we can thus write 

c 

J WP v — (« + l) E l n K ) 

ft 

ft — ft nytl -c{p-Vt) ±^_Y" adf-Vt) 

r 

o- — ft „»T7W , e c(j.-vt) — J^L yw.cO>-w) 
■"0 WP f-(n+l)t — i n e , 

where F n ' , F^', F(" are surface harmonics. 

These do not satisfy the equation of continuity, but the series of derived 
vectors F m , G m) H m satisfy both the equation of continuity and the equations of 
light. 

One of the best forms for these spherical harmonics is to assume one of the 
axes in each of them parallel to one of the coordinate axes. We can then write 



dV 

ft — ft n ay -» c(p-Kt) 

^—°*9 dy S 
dV 

TJ — ft n u Y ~n c(o-Vt) 
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In this case the vector represented by the components F , G , ff is perpen- 
dicular to the surface V_ n = constant. 

Now the components of this vector do not satisfy the equation of continuity, 
and so, although it may represent the vector potential on the electro-magnetic 
theory, it cannot represent either the electric current or magnetic force, and is 
almost without direct use on the elastic solid theory. Hence the derived vectors 
are of more general use. These are 



^i-cO„_ lP jy dz * dy je 
tn-co^p |z dx x dz je 

^i- co »-ip \ x dy y dx \ e 






dJ 



dz 

Where 7 _ f^ + ^L _i_ ^_ _ n tt rf S /r n » E c( P -r« » 

t/o_ dx ^ dy * dz ~ J V ~ n dp *°"P S > 

J =-ncG n _ l9 n - 1 V_ n e^- yt \ 
F 3 = — c*F lt etc. 

The remaining derived vectors are simply repetitions of these with only 
change in the constants. 

These derived vectors can represent the vector potential, the electric 
displacement, the electric current, the magnetic displacement, or the magnetic 
force in the electro-magnetic theory, and the displacement, the velocity of the 
particles or the rotational displacement or velocity of rotation on the elastic solid 
theory. 

The most interesting case is that of symmetry around an axis, say X. In 
this case y _ _<^_ J__ J_ (—1)" 

~ n ~ dx n p~ p n 1.2.3... n^- 1 ' 
where Q n is a zonal surface harmonic. 

If a is the angle measured around the axis of x from y and 6 that made 
with the same axis, and we put p sin = r, 
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we can then write, when m is even 



B m =VG> m + Hl, 



M m+1 =N m+1 r 



1 cOfm-l 



m-2 ^ — r da . ' 



P m =-A a P ro _ 2 = + - 
H m — — A 2 jK m _ 2 — 

_ df. dB,_ ,. r , -ZV m -i_ 1 f <P3C-i , <P3f— i 1 <2ifm-i ] 
m+1— dr "•" da; ~~ ^ iV -»-i~r r * ~ r j tf t tf r dr J" 

The equation of continuity becomes 

d(F m r) d(R m r) _ Q ^ 
dx dr 

These can be expressed in terms of p and $ , as follows : 

Let @ m , P ro and N m be the components in the direction $, p and a respect- 
ively. Then we have 

d . . d cos sin d d 
-= sin 6—, 



dr dp p dpi 

d . d , sin 8 d 

— - = cos 6 -j- H -j- > 

cte a/) /> dp. 

® m =—F m sin6 + R m cos6, 
P m =F m cos6 + R m sm6, 

1 dJf,,.! 



P = — 



p sin dp 
jO 8 d// 



© m+ i=0, 

„ 1 f JJf w _ t BJD , gd»Jf,._ 1 ] 

m+1— josin^l d^ "*" /> 8 d/* 8 J' 

The equation of motion of light becomes in this case 

1 d» Jf m _ x _ d 8 Jf m _ t sin 8 g <f Jf,,.! 
F 8 dtf 8 dp* "*" jo 8 d// 8 ' 
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whence , r AT . a 1 <?Jf m _i 

Now, in this case of symmetry we readily find 

dV 

rr _ ft n uy -n e<j>-Yt) 

dV 
■«o— Wp dr e > 

and for convenience put Z7= C^p'V ' - ™, 

we readily find xt—at . ,. dZ7 f . a^-n a 

17 if 1 = iv 1 r = p sin -r- j sin — j-^- — cos 



dV- n ] 
dr I 



Putting ni _ dQ n _ x 

and, omitting the constant term, we have 

Hence we have, by replacing n — 1 by n, and changing the sign of all the terms, 

M 1 — — G n sin 2 6 Q' n e c{ <- n \ 

From which we find, writing P for the radial component, 
©i=0, 

P 1= 0, 

Ni= — -^sin0g; e c <"-^, 

We can also go backward and write 



p A 

M -i- -#-• 

Hence, as I have remarked before, each of the vectors indicated by these 
components can be made equal to some one of the vectors used in the theory of 
light. Thus 
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Electro-Magnetic Theory. 



Elastic Solid Theory. 



Linear displacement of particle. 
Rotation of particle. 
Linear velocity of particle. 
Angular velocity of particle, etc. 



Vector potential. 

Electric current or electro-motive force. 
Magnetic force, or induction, or mag- 
neto-motive force.* 
Rate of change of current, etc. 
Rate of change of magnetic-force, etc. 

Should we add another term with — i in the place of + i, we shall have 

the following real forms, since 

O n = A n — iB n , 

C' n = A n +iB n , 
c = a — ib, 
d '=■ a-\- ib. 
Ni __ 2smdQ> n ga (? _ n) j ^ cQs b ^ _ _. ^ _ ^ g . n ft (p — -pif) 1 1 



P,= - 



2n(n + l) 
_____ 



Q n e a «-™{A n cosb( 9 -Vt)-B n smb(p-Vt)\, 



N, 



2smdg n (a — ibf 
P 



And the other terms are readily obtained from these as follows : 

N ~i = lta-lf e<lU " n) * [(a ' ~ J2) An + 2abBn] C ° S b ( P ~ Vt) 

— [(a 2 — W) B n — 2abA n ] sin b (p — Vt) } 

e ao>-rt) j ^ a i_ h ^ An _ 2abB n ~] cos b (p - TO) 

- [(a 2 - & 2 )2?„+ 2a___] sin J(p - TC)}. 

As these quantities constantly return again to the same form with only a 
change in the constants, there will be only two cases, when the vector potential, 
or its analogy in the elastic solid theory, is made equal to the even or odd vectors. 
Let us take the first case, making the components of the vector potential equal 
The magnetic induction will then be N x and the components of the 



to © and P 
electric current 



6 > and P « 



4n/x 



4nfx 



* As Mr. Bosanquet has shown a disposition to claim this term, I may say that the idea expressed by 
it is of common occurrence in the works of Faraday, and I myself have used the term and given it 
mathematical expression in the American Journal of Mathematics. Mr. Bosanquet 's article appeared in 
the Phil. Mag. for March, 1883. 
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Hence, in this case, the magnetic induction is in circles around the axis, and the 
electric currents in planes passing through the axis. At a great distance from 
the origin, the electric currents are on the sphere, but the normal component 
must still exist in order that the equation of continuity may be satisfied. The 
case of plane waves is the only one where the normal component entirely 
vanishes. 

"When we make the vector potential equal to N lt the electric currents are 
in circles, and magnetic induction in planes passing through the axis. 

Let us consider the first case. We shall find it much more convenient to 
continue with the imaginary form, remembering that a similar term with i 
changed in sign is always to be added except where it has already been done. 
We shall then have 

Components of the Vector Potential 

^-^^Kn + ^O^ + nC^l^-^ 

_ n(n+l)C n (p _ n) 

Magnetic Induction 

#!=— -^sinfl^o-™. 

Components of the Electric Current 

6> a csmdQ' n i( n , 1 \(Q ,U K (7 ,,!««('-") 

P, _ n(n+l)q, n (p _ TO 



Components of the Electric Displacement 
K d6 __ VKsindQi 



4n dt 47r{2n+l)p 



{(n+ !)<?„_! + nC n+1 \e e <>-™ 



_K_dI\_ VKn(n + l)C n t) 

4n dt ~ ±kc(? Vn 

where (i is the magnetic permeability and K the inductive capacity, and we have 
the well-known relation of Maxwell, 

Y % — — — 
ftK 
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Should a perfectly conducting sphere exist, the electric density on its surface 

will be _ ^ j£_ d_ dM_, _„ _K_d?o 

~ 4jTj0 2 dt dfi 4x dt 

This is the same as the radial component of the electric displacement at this 
point. 

AKBITEAET DI8TOBBANOE. 

In the equations of p. 368 make a = and n = 1. We then have 



tf 2 _tf |l— ^-jy}- 



The magnetic induction is then 



sin ' 



and the components of the electric displacement, 



p, — , t r-o. out. (/ ~ tblf) _ rt) 



i VK cos d 

2nbp* 

Let now a sphere of radius R be circumscribed about the origin and an 
arbitrary uniform displacement take place in the interior of this sphere of a 



value eaual to iVK 

H j£i e -a>(B~vt) — q tvJ ^ £ -ib(R-rt) 



&r6I? l: I 1 bR\- 



Taking this value for the displacement inside the sphere and the previous values 
for the outside, the equation of continuity is satisfied for the electric displace- 
ment and for the magnetic induction. If the sphere is very small indeed, we have 

iVKG 

~ 27Tb 2 R 3 ' 

Whence we have, on substituting the value of (7 from this equation in the others 
and replacing j nR z by dv , the magnetic induction 

M — 2VK P <7 aV ' 

Vol. VI. 
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and the electric displacement, 

&= _ 20+Q |^ sin 6E -iH,-n) dv 

p _ Q ^cos 6e- ib *- rt) dv. 
(7 4^ 

These equations give the complete solution of the problem of finding the 
disturbance at any point due to any arbitrary electrostatic disturbance throughout 
the medium. And it is to be noted that these equations are rigidly exact for all 
distances from the disturbance, and only have to be integrated to give the effect 
of any disturbance. 

Had the disturbance been magnetic, we would have had for the magnetic 

induction, &l= ibpi {20()+ Gi)e -m,-vt ) 

3p 

p„__ 2 cosfl p. £ _ ib(e _ rt) 

P 2 1 

And for the electric displacement, 

N ,__ iVKb sin 6 G e _ i6(p _TO 
Anp l 

Hence, taking a small sphere as before, the magnetic induction within it must be 

* £ ~ r* V bn) 6 

Whence, as before „ 36 ^ r .. , 

in this case for the magnetic induction 

X»ib> S md2C +C i _ v 

~ 8np C„ £ a ' 

SX"ib cos Q iH9 _ rt)fl 
r - 4^~~ G £ aV ' 

and for the electric displacement 

WVKX» S md Q _ vt) 

32h*/> Co 

These equations give the complete solution of the disturbance throughout the 
medium due to an arbitrary magnetic disturbance at any point. 

Although the disturbance is harmonic, yet we know by Fourier's theorem 
that any disturbance can be represented by a series of harmonic terms with the 
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proper coefficients, and, indeed, we can replace the harmonic term by any function 

of p — Vt or o + Vt. Should the disturbance not be parallel to the axis of X, we 

merely have to divide up the disturbance into its components and compute the 

effect of each and then add the components of the computed disturbance. In 

this way we readily find the effect of a general electric or magnetic disturbance. 

Let the components of the arbitrary electric displacement be 

X'e ibrt , X"e ibn , 

Y's ibn , and of the magnetic induction Y"e ibn , 

Z's ibn , Z"e ibvt , 

where, in general, we must replace X', F, Z', and X", Y", Z" by a quantity of 

the complex form, and add another quantity with — i in place of + i. Putting 

D + iE for X', and adding the other term, we would have the real form 

(Z> + iE)e lbn + (D — iE)e~ ibvt = 2 {D cosbVt — E sin bVt\, 
which expresses the disturbance in any phase. But this is only necessary when 
we descend to actual calculation. The effect of this general disturbance is then 
found to be : The electric displacement 

F=z 8^7{ X '( 2 ^+ G W~ 3S , C>x-^^lZ"y-Y»ziy< b «-™dv, 

S'= ^ 9 {z\2C + G i )^-SS'C i z- S ^^[Y''x-X''y]y ib ^^dv. 

The magnetic induction, 
F"= ^ S {X"(2C + <7 2 )p 2 - SS"C,x + 1 ^P[Z'y- F«] }*-•&-*»*,, 

G "= 8^/{ F "( 2 ^+ <W~ SiSPCto + 1 -^[X'z-Z'x]y^-^dv, 

where I have written S' = X'x + Y'y -\-Z'z, 

S" = X"x + Y"y + Z"z. 

For a general expression the term 6 - «&>-»'*> must be replaced by any function 
of o — Vt or o + Vt as before. 

Before integration, one must of course substitute (x — x'), (y — y 1 ) and 
(z — z 1 ) for the x, y and z of the formula. It is evident that the components of 
the electric displacement can be replaced by the true convection displacement of 
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electricity as carried along by the actual motion of the medium, and the disturb- 
ance due to a moving magnet can be calculated in a similar manner. 

To obtain an idea of the relative magnitude of the quantities which enter 

into these expressions, I may remark that the value of b in wave lengths is —r~ • 

Hence, for any ordinary calculations with respect to light, the distance need be 
only a few inches, or, indeed, one inch, to cause the values of C lf (7 2 , etc. to 
become constant and equal to G . But if we are treating of longer waves we 
must retain such terms. 

I have already given the proper expansion into series of the quantities here 
involved, but it will be better to put them in a reduced form, 



c = — ib ; p=p */l -\- 2s ; cps = q . 

^ = £{fll<ri + { 0,( P ) + & %,) + eta} 
^±^'=f{i(^. + tf,W) + etc. 



P 



+ 



[(2n + 2)G n (p) + G n+i (p)-] + etc. j 



(2n + 3)1.2.3. 

Thus one has p 2 = IP — 2 (aa/ + ytf + «0 -f A 

where R*= x 2 +£/ 2 + £ and r* = oP -f y>*+ g' 2 . 

Hence one can write „ xri -\- ytf -\- zd r 2 

p-M; s- — ^ +2E 2 ' 

xx? + yi/ + zz 1 , r 2 

When R is great and r is small, the last term will vanish and leave a very conve- 
nient form for many cases. 

But it is not necessary to always perform the calculation at a finite distance, 
for we have the following theorem which I believe to be new. 

Theorem. Supposing the sources of light to be continuous and knowing the 
light over a sphere at an infinite distance, it can be determined for all space by 
the following process : 

Let F, G and H be the components of any of the vectors such as the 
electric displacement, etc., and suppose them known over the infinite sphere. 
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Express them in term's of series of surface harmonics, thus : 

F=^- {E Y ' + E 1 Y^ + E i Y i '+ etc.} 

G = G » t0{p - n) \e q y » + E X Y{' + E,Y» + etc.} 

H=Q^^{E Y '" + E 1 Yr + E i Y i '"+ etc.} 

At any other point of space the equations must satisfy the equations of 
light motion and also the equation of continuity. Referring to the equations of 
pp. 363 and 364, we see that if we multiply each term by the corresponding 
quantity G n , the equations of light will be satisfied, and if the surface harmonics 
are of the form 

the equation of continuity is also satisfied. Hence the components of the vector 
at any point of space are 

F- ^p-{C Q E Y '+ O&Yf + O&Y,' + etc.} 
Q _^p.\c i E^ , + C x E 1 T l H +C t E t Yi'+ etc.} 

H= e -^p-{C E Y,!"+C 1 E 1 Y 1 "'+ C % E % Y1"+ etc.} 

The best point for the origin will be somewhere near the center of gravity of 
the illuminated body. 

That the light is perfectly determined in this way for all points outside a 
sphere around the origin which does not cut the illuminated body is evident from 
the fact that we might reverse the motions so as to make the light return along 
its previous path. 

These expansions and this theorem entirely change the ideas of those who 
have only been in the habit of regarding light from the point of view of rays of 
light. For we here see that light coming from any source may be replaced by 
light from another source, so that the true source might be entirely invisible, and 
we " see " only the false source. This, however, could only happen for a point 
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outside a sphere drawn around the false point and the real point, and we could 
always detect the deception by making a complete spherical journey around the 
real point and inside the false point. 

As an illustration, let us compute the effect of a circular electric current 
which is caused to vibrate back and forth according to a simple harmonic function, 
so as to make the displacement C e~ ert . In this case we must make 

Q = X' = X"=Y"=Z"; Y'dv = —O r sin ada; Z'= C r cos ada, 
where r is the radius of the circle. 

From the symmetry around the axis of x we have F'=0, and can write for 
the component of the electric displacement around a; at a distance R sin from 
that axis, _ _ r Wtr cn (*> f 2O + Q „.„ ZRrQ> sin 6 



-I < — — — -sin a \ cos*a|-e'" , aa, 



Whence ■£- (^kt3 e ^ = - Rr s i n C os a -i ^ [4 G x + <7 8 ] , 



8tt 

p 2 = iE 2 — 2i2r sin sin a + r 2 . 
Integrating the first term by parts, taking sin ada for one part we have, 
since the first part disappears, 

jf fc *y^ liBofc =jp£(!«±s,.) cos**,. 

d dp d Rr sin cos a d 

da da dp p dp 

1 - 1 — -V p ) = — i2r sin0 cosa-^-^ 
p J 5 p* 

Whence we have 

N'= s I —V- cosrada. 

8n *n p- 

In the expansion, put 

r 2 
p-=R; q= — cr sin0 sin a + c-^p- 

~ = 5 { W + T ^ (i2) + A ^ + etc -} 

Writing = £andA = -ersin0 f 

we have # = g + h sin a , 

Jo 2 cosa « a — ^{ 2 # T 12 2A 9 « i- 1.2.3.4 2.4.6^ 

4. etc + "("-!)• ••("-2m+l) 1.3... (2m- 1 ) m J 

+ etc - + 12... 2m 2.4...(2m+2) 5 ' * i" etc. j 
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So that the problem is completely solved without any approximation and for all 
distances at which the series is convergent. 

At a great distance the expression becomes very simple. 

iV 4tt B e \2^ 2?A arS U ^ (2.4) 8 6 ° T P 

+ GtC - + (2.4...2n)^ + 2) ^* ^ + etC '} 
We recognize this series as a Bessel's function divided by br sin 6 , and so we can 

write , T/ 8ib?rC e(R _ rt) T ,, • m 

iv — e c iX rt> J x (br sin 6) . 

This same series occurs in the expression for the light from a circular orifice, but 
I am not aware that writers on physical optics have recognized this connection 
with Bessel's functions. Prof. Stokes has given the value of the series for large 
values of br sin 0, and the same value is given by writers on Bessel's functions. 
We thus find, writing v = br sin 






3.5.1 1 3.5.7.9.1.3.5 1 



sin (v— ^j\ 1 + : ™ —-.,— "-Tr"^ ~~, + etc. 



2.4 (4t>) 2 2.4.6.8 (4v)< 

/ tt\(3 1 3.5.7.1.3 1 , t 
| + C0S ^__)J__ ____ +etc 

The energy given out per unit of time is, by the first formula, 

2K 1 5 + 35 720 + 31680 1310400 + J 

and by the second for very large values of br 

VrWCl 

27r.gr 

Let us now reduce the series to spherical harmonics. 

Write ., 2nir 

u = cr — — %br = — > 

A 

P = W lT + W sin20 + (2^6 Sill40+etC -} 

This satisfies the differential equation 

ePP 1 dP . 

-3-5- 3 P sin 2 0=0. 

air m aw 

Write also /« = cos and (?„ = ~ . P can be developed in the following series 

as it only contains the even powers of ft ; 

P=A\ + A' 3 Q> 3 + AW* + etc. 
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But 

"Writing o = -,— r -3— we have 

aw* w aw 

9.11 f 3.4 2 _ | 

^ 5_ 5.6 r 5. 9^ 3 15 ^ x dA *y 

A' 1 = etc., 

,, _ (2"-l)(2n + l) f (n-2)(n-l) ,, (»- 4 )( n -3) ,, _.,, 1 

^»— n(n + l) I (2n— 5)(2w — 1)"*— » (2n — 5)(2n— 7) n ~ 4 0J1 — »J* 
These all depend on A\, whose value can be found by developing sin 2 , sin 4 , 
etc., in harmonics. In the following values one must add the second term with 
— i in place of + i before substituting in the formula 

A'i=-§-<W« ,r 

91 



^'■= § {t <«••) + a °m4 



A\ = etc. 
These reduce to Bessel's functions when we obtain the real form. Thus 

Each term of the original series is now of the form A' n Q n sin 0. But by the 
theorem of p. 372 we have only to multiply the terms of this by C n (R) and it 
gives the value for any finite distance outside a sphere around the origin which 
contains the circle. Hence 

N ' = ^&W e c{B - n) \ A 'i°o+ A', q 3 G 3 (B) + A' t Qf s G t (B) + etc. f . 

Dynamical Theoby of Diffraction. 

When a ray of light strikes upon a screen in which there is an opening, a 
disturbance takes place in that opening, whose effect can be calculated by the 
preceding formulae. Maxwell has shown, that in a plane wave, the energy is 
half magnetic and half electrostatic, and that the magnetic and electric displace- 
ments are at right angles to each other and in the same phase. 
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Stokes' solution is based upon the displacement and rate of displacement of 
his elastic medium. But in an elastic wave there is not only displacement but 
rotation also, and the components of this rotation must satisfy the equation of 
continuity. But when a wave is broken up at an orifice, the rotation is left 
discontinuous by Stokes' solution, and hence it cannot be exact. The equation 
of propagation of the rotation is the same as that of the displacement, and the 
two are at right angles to each other, and they are both equally important. 

Hence, on the elastic solid iheorg as well as the electro-magnetic theorg, the true 
solution of diffraction will depend upon the sums of two similar terms. 

When we take account of both the terms, the fundamental properties by 
which Stokes attempted to obtain the direction of the motion of the particles 
vanishes, and that problem is impossible of solution by this means. 

In the general equations of p. 371 let all the disturbances vanish except X' 
and Y", so that the electric disturbance is in the direction of X and the magnetic 
in that of Y. 

In order that the energy coming from the two disturbances may be equal 
we must have V^K / WX' \*_ V*K / 3b 2 Y" \* 

8tt \2VKj~ Stt \ 8x J ' 
or v „_ 4?r 

Y ~Tk x 

The electric displacement at any other point of space will then be 

WY'v 

F=-l^\(iO t +C t )f?-SO, a ?+SO, 9 zle-""-">d s , 

IfiY'v 

H'=^4ril—^O i xz-3G 1 xole- ibf <'- n) ds. 

For showing the peculiarities of the case, polar coordinates are best. Let be 
the angle made by p with the axis of z and <£ the angle around it from the plane 
XY. Let &, 0", 4>', 4>", and P, P' be the components of the electric displace- 
ment and magnetic induction to increase these angles and in the direction of o 
respectively. Then we have for the electric displacement 

^. W*X'v f r, . /nT, * ~\ COS0) ... ..,. , 

B'=: ~ 7 , -COS4>|[l + COS0][l--^]--^| e - 6 "'-^ < fe, 



<D'= — 



illp 

WX'v 



Sxp 
SibX' u 

4:7Tp 2 



sin t {[1 + cos 0][l — JL] - g^} *-»*-*>&, 



P~ "^ly sin 6 cos $ 1 1 — ^-1 e-^-^ds. 
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P" - M ^ v sin sin $1 1 — ~ le- ib{ "- n) ds. 
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And for the magnetic induction, 

*"=w cos *i [1+cose] [ 1 -^-^?} 8 """"'"'*' 

MbX'v 
VKp 

In these, v is the thickness of the disturbed stratum, so that vds = dv. The 
electric displacement within the stratum is X'e ibvt . 

When bp is very large, the disturbance is in the spherical surface, and indeed 
it only requires a fraction of an inch from the origin to be able to omit the 
terms in bp, and also P and P". The equations then become very simple, as 
follows. The electric displacement, 

0' = ¥£l cos $ (1 + cos 6) e- ib «- n) ds , 

q I= O0 ± V g x + COB m 6 -i»(p-F») (fo| 

OTZO 

and the magnetic induction, 

P"= |^ cos*(l + cosd)s- ib( "- vt) ds. 

2 VKp r v ' 

From either of these expressions or the more exact ones we see that the 
distribution of the magnetic induction is exactly the same, but turned around 
the axis of Z 90°, as the electric displacement. As this result would also apply 
to the elastic solid theory, we conclude that diffraction gives no means of deter- 
mining the relation behveen the displacement and plane of polarization. Had we 
only taken the original disturbance as electrical alone, we would have arrived at 
Stokes' result. 

Squaring the coefficients of the time function and adding, we find 

Hence the light is symmetrical around the axis of Z and varies from 4 in the 
positive direction to 1 in the plane XY and in the negative direction. 

It now remains to connect the arbitrary displacement X'v with the intensity 
of the original wave. 
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Let the arbitrary displacement X've iin , with the corresponding magnetic 
quantity, exist throughout the plane XT. From considerations of symmetry, 
the electric displacement throughout space must be parallel to the axis X, and 
so we can write, 

z— constant, 

«= — rcos^>, ds = rd<j>dr = pd<pdp , 
y = — rsin<£>. p 2 =z 2 + ^- 
The general and exact integral is 

J>= saz>vf 1 i+2b* * tf i _ iMp _ ro 

8 \ bp p 2 bp 9 J 

For positive values of z this gives, between the limits p = oo and p = 2, 

2 

But for a negative it is zero. Hence such an arbitrary disturbance produces a 
wave in the positive direction, but none in the negative direction. 

No approximation has been made in obtaining this quantity, and it evidently 
applies to a plane of any size, even infinitessimal, provided the point under 
consideration is infinitely near to it. The displacement near the surface therefore 
differs in phase -| n from the arbitrary disturbance, but is dependent upon its 
value at that particular point. 

Hence we can replace any particular wave surface by a surface of arbitrary 
disturbance whose phase differs \it from that of the wave. Such a surface of 
arbitrary disturbance then produces the same effect at all points of space as the 
original wave. Should the wave be spherical and of short radius, the normal 
component will enter and complicate the result, though the solution can be 
obtained but would evidently be complicated. 

But if the wave have a radius of an inch or even less, the displacement is 
practically perpendicular to the radius, and the solution here given will apply. 

Let I, m, n be the direction cosines of the normal to such a wave, and V, tn! 
and n' the direction cosines of the electric displacement which is represented by 

Then the values of the arbitrary displacement and magnetic induction to substi- 
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tute in the general equations of p. 371, to produce such a wave will evidently be 

Where Z", wi", »" are the direction cosines of the magnetic induction. We also 
have IV + mm' + /m'=0, 

ZZ" + m«z" + «w" = , 
ri"+m'm" + rin"=0. 

Returning again to the equations of p. 377, we see that the polarized light 
is diffracted equally in all directions from a very small orifice and independent 
of its plane of polarization. Furthermore the plane of polarization at any point 
is found by drawing a sphere through that point with its center at the orifice, 
and then drawing a plane through the given point and the point where the inci- 
dent light first cuts the sphere, and cutting the plane of polarization of the inci- 
dent light in a line perpendicular to the incident ray. The intersection of the 
plane and sphere then give the direction of the polarization. 

It is seen that in both these particulars my solution differs from that of Prof. 
Stokes, and the construction is the same whether one takes the electric or mag- 
netic quantities as the direction Of polarization. 

The system of planes for the electric and magnetic quantities form a system 
of orthogonal circles on the sphere. 

The following constructions can also be used for obtaining the direction of 
the electric and magnetic quantities. Draw a sphere around the orifice and draw 
an axis through the sphere in the direction of the incident light. Then rule a 
sheet of paper with a series of lines at right angles to each other. Cut a star 
shaped piece out of the paper with its diameter equal to the circumference of 
the sphere, and having a very large number of points. Place the center of the 
star on the sphere at the end of the axis where the light leaves the sphere, and 
wrap the points around the sphere, the points meeting around the incident ray. 
The marks on the paper then give the required directions. 

We can also construct the curves of polarization by noting that the stereo- 
graphic projection of the lines on a plane is merely a series of straight lines. 
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The equations become very simple at many wave lengths distance from the 
orifice, especially when the orifice is small. If the radius of the original wave 
is large, it is usually sufficient to consider the periodic factor as the only variable. 
In this case we can write 

&=-£f cos 4> [1 + cos B~\ffe- ib ^- n) ds, 






$' = — — s in$[l + cos eiffe-w-^ds, 



P'=0, 

where / is the coefficient of the original vibration, and therefore its square is the 
intensity of the original light. The intensity of the diffracted light is simply 
proportional to the sum of the squares of 0* and <&. This is the expression 
ordinarily used except the term in 0. Thus for a circular orifice we have the 
vibration expressed in Bessel's functions, 

0'=^cos^(l + cos0)^t^- ) e--'^', 

tf= —g; sin<Kl + cos6) Mb : 8i » d) e-™-™ 
&K sin a 

It is impossible to pass from these expressions to the case of a plane wave 

since they are only for the case of a great distance from a small orifice. 



